Spin dynamics in InAs-nanowire quantum-dots coupled to a transmission line 
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We study theoretically electron spins in nanowire quantum dots placed inside a transmission line 
resonator. Because of the spin-orbit interaction, the spins couple to the electric component of the 
resonator electromagnetic field and enable coherent manipulation, storage, and read-out of quantum 
information in an all-electrical fashion. Coupling between distant quantum-dot spins, in one and the 
same or different nanowires, can be efficiently performed via the resonator mode either in real time or 
through virtual processes. For the latter case we derive an effective spin-entangling interaction and 
suggest means to turn it on and off. We consider both transverse and longitudinal types of nanowire 
quantum-dots and compare their manipulation timescales against the spin relaxation times. For 
this, we evaluate the rates for spin relaxation induced by the nanowire vibrations (phonons) and 
show that, as a result of phonon confinement in the nanowire, this rate is a strongly varying function 
of the spin operation frequency and thus can be drastically reduced compared to lateral quantum 
dots in GaAs. Our scheme is a step forward to the formation of hybrid structures where qubits of 
different nature can be integrated in a single device. 

PACS numbers: 73.63.Kv,73.63.Nm,72.25.Rb 



I. INTRODUCTION 

Over the last decade, the spin of individual electrons 
in semiconductor nanostructures has been intensively 
studied in relation to spin-based quantum computing 
schemes Attaining an almost full control over the 
spin of individual electrons in QDs opens the possibility 
to study single spin dynamics in a solid state environment 
in the presence of relaxation and decoherence. Although 
lateral QDs have been most successfully used until now 
to demonstrate spin coherence and usability for quan- 
tum computing*^ novel quantum systems have emerged 
in recent years, providing a number of new ways to im- 
plement the basic ideas of quantum computing, 4 Among 
such systems are the QDs formed inside semiconductor 
nanowires. 5 '^ 

Rapid progress in GaAs nanostructures started once 
few-electron QDs became available (for a review, see e.g., 
Ref0), which opened the door to control the number of 
electrons in a single QD down to one in vertical^ and 
lateral 9 dots, as well as in double QDs . 10 ! 11 ) 12 Further 
important experimental progress came with the advent 
of charge sensors which, quite remarkably, enabled the 
measurement of the relaxation time of one single spinJ-S 
The longest spin relaxation times in single GaAs QDs ex- 
tend up to several seconds^ - and were measured in lateral 
dots at relatively small magnetic fields (B ~ IT). 

The spin decoherence time in GaAs was measured also 
in double QDs by studying the hyperfine-induced mix- 
ing of singlet and triplet statesi 15 i 16 In the same set-up, 
a universal entanglement operation was implemented^ 
enabling a square-root-of-swap operation 1 , between two 
spin-1/2 qubits on a time scale of 180 ps. Resonant 
and coherent manipulation of a single spin-1/2 has re- 
cently been implemented in a GaAs double QD, mak- 
ing use of electron spin resonance (ESR)^^ as well as 
electric dipole induced spin resonance (EDSR) 19 i 20 tech- 
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FIG. 1: Schematics of the two configurations considered in 
this work. A) Large-diameter InAs nanowire (pink) posi- 
tioned inside and parallel to the transmission line (blue). The 
disk-shaped quantum dots (QD) are located in the nanowire 
and are formed by two InP-boundaries (brown). Each QD 
contains only one electron with spin 1/2 (green). B) Two 
small-diameter InAs nanowires (pink) positioned perpendicu- 
larly to the transmission line (blue). The elongated QDs are 
oriented along the nanowire with one electron (green) in each 
dot. The QD confinement can be achieved by barrier materi- 
als (as shown in brown) or by external gates (not shown). 
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spin manipulation in double dots was also recently 
demonstrated^ 

The use of different semiconductors, other than GaAs, 
has since long been pursuit with the goal to create nanos- 
tructures with novel properties. Particular examples are 
InAs and InP nanowires, where both gate defined and 
'barrier' defined QDs could be fabricatediS2i22i2ijS& The 
advantage of these materials is that both optical and 
transport measurements can be carried out on the same 
type of structure. The number of electrons can equally 
well be controlled down to one electron per dotj2£ which 
shows that QDs created in nanowires can serve as alter- 
native candidates for spin-qubits. 

One particular difference between GaAs and InAs 
semiconductors is the strength of the spin-orbit interac- 
tion (SOI), which is much larger for the latter material. 
This fact, however, is a double-edge sword; on one hand 
it opens up the possibility to efficiently manipulate the 
electron spin with electric fields only , 20 ' 26 i 27 ' 28 i 29 while on 
the other hand it implies stronger coupling of the spin to 
charge environments, like phonons, particle-hole excita- 
tions, gate voltage fluctuation, etc. However, due to the 
quasi-lD structure of the nanowires, the spin relaxation 
times due to phonons and SOI turn out to be longer 
than one might expect from QDs created in InAs bulk 
material. Indeed, the time scales obtained in this work 
are on the order of microseconds to milliseconds for suf- 
ficiently large Zeeman splittings. At the same time, the 
relaxation rate exhibits peaks as a function of a static ap- 
plied magnetic field due to the quantization of the phonon 
spectrum. The long relaxation time and the presence of 
strong SOI permits then an efficient control of coherent 
spin states by making use of EDSR^^^^^ 

One of the main ingredients in the spin-qubit scheme 1 - 
is the electrical control of two-qubit gates to gener- 
ate entanglement. While the original proposal involved 
only local interactions between neighboring spins, it is 
desirable to couple spins directly over large distances, 
since this produces a better threshold for fault tolerant 
quantum computation^ A solution to this problem was 
first proposed in Ref. [3l| and involves optical cavities 
whose photon modes mediate interaction between dis- 
tant spins. The coupling of the spin to optical cavities 
in semiconductors was also the subject of some recent 
experiments! 32 ' 33 

Very recently, ID electromagnetic cavities (or trans- 
mission lines) were shown to be very suitable for reach- 
ing the strong coupling regime between superconducting 
qubits and photons i 34 ' 35 ' 36 Theoretical extension to QDs 
were proposed subsequently, including charge and spin 
qubits j 37 ' 38 The direct coupling of the spin to the cav- 
ity modes via the magnetic dipole transitions is usually 
weak and one has to use electric dipole transitions to- 
gether with correlations between spin and charge degrees 
of freedom in order to obtain a strong effective coupling. 
This can be achieved in several ways, e.g. by making use 
of the Pauli exclusion principle and Coulomb repulsion, 38 
or of Raman transitions^ 7 - 



Here we propose another mechanism to achieve long- 
distance coupling between spins inside a cavity, namely 
via SOI which leads to an effective coupling of spin to 
the electric field component of the cavity photon, and 
thus eventually to a coupling between distant spins me- 
diated by this photon. In order to reach a sizable cou- 
pling strength, it is desirable to use nanostructures with 
large SOI such as InAs QDs. Two such proposed con- 
figurations, which define the two model systems to be 
studied in this paper, are sketched in Figs. Q]A and B. 
They consist of nanowire QDs embedded in a transmis- 
sion line. In particular, in Fig. [T]A_ a nanowire positioned 
parallel to the transmission line axis is shown. In this 
case, the QDs are realized by confining the electrons in 
the longitudinal direction (i.e. along the nanowire axis) 
much stronger than in the transverse one. This corre- 
sponds to a nanowire with a large diameter, on the order 
of 80 — 100 nm. Such longitudinal confinement can be 
achieved by applying metallic gates or by using other 
materials as barriers (InP for example, which is depicted 
in Fig. [T]in brown) which have a larger band gap than the 
host material such as e.g. InAs*^ 3 - In Fig [TJ3 a small- 
diameter (D < 40 nm) InAs nanowire is shown, being 
positioned perpendicularly to the transmission line and 
containg QDs that are elongated along the nanowire. 
That means that in this case we assume that the elec- 
tronic confinement along the nanowire is much weaker 
than in the transverse direction. Then, to a very good 
approximation, the electrons can be considered as be- 
having one-dimensionally, which will allow us to treat 
the SOI exactly, while this is not possible for the config- 
uration Fig. [T]A However, in order to prevent a current 
flow, the nanowire and the transmission line need to be 
separated by some insulating coating material, obtained, 
for example, by atomic layer deposition. 

The goal of our work is now to analyze these config- 
urations in detail and, in the first part of the paper, to 
derive an effective spin-spin coupling Hamiltonian. In 
the second part, we study the spin decay in this system, 
induced by phonons and SOI, and calculate explicitly the 
spin relaxation and dccohercncc times due to this mecha- 
nism. We will show that these times are much longer than 
the switching times needed to manipulate and couple the 
spins coherently. Thus, our findings provide theoretical 
evidence that nanowire QDs embedded into transmission 
lines are promising candidates for spin-qubits with tun- 
able long-range coupling. This scheme also opens the 
door to hybrid configurations where qubits of different 
nature (e.g. superconducting and spin qubits) can be 
coupled via the transmission line. 

The paper is organized as follows. In Sec. II we in- 
troduce the model for our system, namely single-electron 
QDs and cavity and specify the model Hamiltonian. In 
Sec. Ill we derive first the effective spin-photon Hamilto- 
nian for a single spin in the cavity for a general SOI. Here 
we derive also the general effective spin-spin coupling in- 
duced by the SOI and the cavity photon modes. In Sec. 
IV we investigate the case of a QD strongly confined in 
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the longitudinal direction. Then, in Sec. V we analyze 
the opposite case, when the electron is strongly confined 
in the transverse direction of the nanowire. In Sec. VI 
we provide some numerical estimates for the strengths of 
the spin-photon and spin-spin couplings for both cases. 
Then, in Sec. VII we give a brief description of the ma- 
nipulation of the spins by electric fields. In Sec. VIII 
we study the spin decay and provide a detailed descrip- 
tion of the relaxation of the spin via SOI and acoustic 
phonons. Some technical details of the phonon analysis 
are deferred to App. A. Finally, conclusions are given in 
Sec. IX. 



II. MODEL HAMILTONIAN 

The Hamiltonian of the system composed of the single- 
electron QD and the cavity reads 

H= ^ + V(r) + ^g f x B B-a + Hso + H e _ J +H^ (1) 

where the first two terms represent the bare orbital part 
of the Hamiltonian, m* is the effective mass of the elec- 
tron, g is the g-factor of the electron in the material, and 
V(r) is the confinement potential, both in the longitudi- 
nal and transverse directions. We can obtain an effective 
Hamiltonian H e g by averaging over the ground-state |0) 
in the longitudinal or in the transverse directions depend- 
ing on which case in Fig. [T] is considered. Then, for the 
system in Fig. Q}\(B) we obtain an effective 2D (ID) 
Hamiltonian. 

The third term stands for the Zeeman interaction, 
while the fourth term in Eq. ([T|) represents the SOI. For 
wurtzite InAs nanowires grown along the c-axis, with the 
longitudinal confinement much stronger than the trans- 
verse one (see Fig. [T]4) the SOI takes the form of a 
Rashba type^ Hso = H so = a (P x c ) ' a i which, when 
written in components, becomes 

H so = a(p x a y - p y <T x ) . (2) 

We mention that our present study is quite general and 
can be easily adapted to other types of SOIs (such as 
Dresselhaus type). In the opposite case, when the trans- 
verse confinement is much stronger than the longitudinal 
one (see Fig. [TJ3), the SOI Hamiltonian Hso takes the 
form Hso = H l so = (k ■ c)(tj ■ a) which, when written in 
components, becomes 

H SO = VPxfJr,, (3) 

with r\ = (rjx, T)y,f]z) being a vector of coupling constants 
and a v being the spin component along r/. 5 

The sixth term represents the interaction between the 
photons in the cavity, labeled 7, and the electron in the 
QD. This term is given by 

ff e _ 7 = eE(z) ■ r. (4) 



The electric field E(z) acting on the electron is E(z) — 
e y V(z)/d, with e v being the unit vector along y, V(z) 
represents the fluctuating potential within the transmis- 
sion line and d is the distance between the transmission 
line and the center conductor. The voltage fluctuation 
V(z) has the following form2£: 

00 

V(z) = J2v P sm(^)[a p + al], (5) 
P =i 

where V p — ^Tiuip/Lc, a P {a v ) are the creation (anni- 
hilation) operators for the excitations(photons), c the 
capacitance per unit length, L the legth of the res- 
onator, and ujp the eigenmodes of the resonator. The 
last term in the Hamiltonian represents the free photons 

From Eq. {T]) we see that there exists an infinite num- 
ber of frequencies in the transmission line, implying a 
coupling of the electron charge to an infinite number of 
modes. However, from all these modes, the relevant ones 
are those close to resonance with the Zeeman splitting of 
the spin. In the following we disregard all other modes 
from the problem and we assume also that the QD is in 
the center of the transmission line, so that the interaction 
between the electron charge and the photons becomes 
maximal. Having now defined all the ingredients, we can 
proceed to study the dynamics of the system. 

III. GENERAL SPIN-PHOTON DYNAMICS 
A. Spin-photon interaction 

In the following we derive an effective spin-photon 
Hamiltonian, assuming for both cases in Fig. Q] a SOI of 
arbitrary strength (to be restricted later on) . In the case 
of time-reversal symmetry, the ground state of the dot 
{Hd = H + Hso + Hz) is two- fold degenerate (Kramers 
doublet), while this degeneracy is lifted in the presence 
of a magnetic field. If the magnetic field is such that the 
doublet splitting and also the electron-photon coupling 
strength are smaller than the level spacing of the QD, 
we can restrict our considerations to the dynamics of the 
lowest doublet only. We label this doublet by {| -ft), | JJ.)}, 
which is now different from the 'true' electron spin. 

We can connect formally the states in the presence of 
the SOI to the ones in the absence of the SOI with the 
help of a unitary transformation or Schrieffer- Wolff (SW) 
transformation 

K) =e- s \n)\a), (6) 

where the states \n) are the eigenstates of the Hamilto- 
nian Hq (Ho\n) = E®\n)), \n T ) are the Kramers dou- 
blets with SOI, \a) = I t) I) are the bare spin states, and 
S = -S^. Also, the relation H d \n T ) = E^ T \n T ) holds 
from our definition of the transformed state. For nota- 
tional convenience we denote the lowest Kramers doublet 
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as |O r ). This is written simply as |O r ) = |r), with the 
identification |t) = {| f|-), | if)}. The above transforma- 
tion can be performed on the level of the Hamiltonian, 
implying diagonalization of the Hamiltonian Hd in the 
basis of the 'bare' Hamiltonian Hg 

H = e- s He s . (7) 

The advantage of transforming the Hamiltonian Hd so 
that it becomes diagonal in the basis of the bare Hamilto- 
nian Hq is now obvious. Within this transformation one 
can in principle proceed to calculate the effect of SOI to 
arbitrary order in perturbation theory, together with the 
SOI induced spin-photon coupling. We can now derive an 
effective spin-photon Hamiltonian within the lowest dou- 
blet |r) by averaging H over the orbital ground state |0). 
This leaves us with the following effective spin-photon 
Hamiltonian i? s _ 7 = (0|H|0) given by 

H s - 1 = -gn B B e ff<j z + M. 1 ■ <r(a t + a) + huja^a, (8) 
where 

l -gix B B eB G z = l$i\er s H d e s \$) (9) 

stands for the renormalized magnetic field and 

M 7 • a = ^(0|e- V|0>. (10) 

stands for the spin-photon coupling. We mention that in 
order to have a finite coupling of the spin <x to the pho- 
tons, the vector A^ 7 must contain some time-reversal 
breaking parameter, such as the external magnetic field 
B. In the absence of the magnetic field there is no 
coupling between the lowest doublet and the photons 
(At 7 = 0) to all orders in SOI. 
We now define the spin-photon coupling strength v = 

<J (M.j) 2 + (M.%/) 2 and the detuning of the qubit from 

the cavity by A = E e z — hui, where E^P = g[iBB e g. 
Close to the resonance between the qubit and the cavity 
mode (A <C E^f , Too) one can simplify Eq. ((SJ by using 
the so called rotating wave approximation (RWA) ^ This 
implies to switch first to the interaction picture, so that 
the operators a(a)) and er^, where er^ = a x ^fia y become 
time-dependent 

a T (t) = a T (0)e^ fft (11) 
a(t) = a{0)e~ wt (12) 
<Jz(t) = tr,(0). (13) 

where ui z = E z /%■ Then, we neglect the terms in 
the time-dependent resulting Hamiltonian which oscil- 
late fast on the time scale %/A. This means neglect- 
ing counter-rotating terms such as aV + ~ e 1 '"*^'', 
aa- ~ e -i(f»f+<*)t^ a \ Uz ^ ^ and a(Jz ^ e -iut^ which 



average to zero for large times. Within this approxima- 
tion the Hamiltonian in Eq. ([5]) within the interaction 
picture becomes static and of the form 

Hs- 7 = -j9l l BB e ff<j z + v{a x a- + a+a) + hiva^a. (14) 

As expected, the above expression agrees with the 
Jaynes-Cummings Hamiltonian 4^ 

B. Effective spin-spin interaction 

We now investigate the case of two QDs in the cavity 
in the limit of finite detunings Ai ^. The Hamiltonian 

(2) 

Hg' corresponding to the two spins in the cavity can be 
found by just extending Eq. (|14[) to two spins 

(15) 

For Vi/Ai < 1 (i = 1, 2), the spin-photon interaction can 
be treated within the second order perturbation theory in 
Vi. We use again the SW transformation, similar to the 
previous section. Here, this implies finding an operator 
T so that 

= e r ifi V T ( 16 ) 
is diagonal in the basis of the spin-photon Hamiltonian 

(21 

without spin-photon interaction (the Hamiltonian H\' 
with v\ 2 = 0). Within first order in spin-photon cou- 
plings V\ 2) the transformation operator T reads 

T = E ^««-«vi), a?) 

i=l,2 1 

under the assumption that the condition fj/Aj < 1, 
i = 1,2, is satisfied for both dots. The transformed 

— (2) 

Hamiltonian H:_^, becomes 

+ {Ef z+ ±).i + {Ef z+ ±y z 

(-^ + ( cr + (T - + oW-), (18) 

where — g^BB l e g . We can obtain a pure spin Hamil- 
tonian by neglecting the fluctuations of the photon num- 
ber a^a — > (a^a) = n, with n the average number of 
photons in the lowest cavity mode. The resulting Hamil- 

— (2) 

tonian H s = flg_ 7 | to-»n reads 

H s = E\g\ + E\a\ + J{aX(j 2 _ + cr+ffi), (19) 

where 

& z = E* + 2(n+~\£, t = l,2, (20) 
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J = uiu 2 



1 



(21) 



In Eq. (|20| we see that the effective Zeeman splitting 
E l z is quite different from the bare one, Eiz- Besides 
the SOI renormalization of the Zeeman splitting, there is 
also a contribution from the spin-photon coupling, which 
consists of the Lamb shift (the term independent of the 
average photon number n) and the ac Stark shift (the 
term proportional to the average photon number n). 

The expression Eq. (fl"9|) is one of our main results: in 
the presence of SOI and cavity modes one can achieve an 
effective spin-spin coupling with the exchange coupling J 
between two spins that are spatially well-separated. In- 
deed, this interaction can act over the entire length of 
the cavity, which can be as large as a few millimeters. 
Also, the spin-spin interaction is of XY-type (transverse 
spin-spin coupling), which together with single spin ro- 
tations have been shown to be universal for quantum 
computing ] 31 ' 41 We mention that in order to obtain a 
maximal effect, one should be able to tune the two qubits 



into resonance, so that E\ 



El 31 



IV. STRONG LONGITUDINAL 
CONFINEMENT 

So far we have taken the SOI into account exactly, 
regardless of the system under consideration, but under 
the assumption that the lowest Kramers doublet is well 
separated from the higher states compared to Zeeman 
energy and electron-photon coupling. We analyze here 
the spin-photon coupling for the case shown in Fig. [TJ\. 
As stated in Section II, in this case we can derive an 
effective transverse Hamiltonian H e g = H t = (0i\H\0i), 
where \0i) stands for the ground-state wave-function in 
the longitudinal direction z. The effective Hamiltonian 
H t reads 



Hi 



pI 



V 



2m* 



v - + V{x,y) + H z + H t so + H e ^ + H 1 , (22) 



with V(x,y) = (0i\V(r)\0i), while all the other terms 
stay the same since they do not act in the z-direction. 
In the above expression we disregarded the term 
(0;|(p^/2m*)|0;), as it gives a constant shift of the levels. 

We can start in principle to derive the spin-photon 
interaction from the effective Hamiltonian H t by mak- 
ing use of the transformation ([5]). However, this cannot 
be done exactly and we have to proceed in perturba- 
tion theory. In order to give some numerical estimates 
for the strength of the coupling u, we assume the limit 
of weak SOI, quantified by the condition R/Xso < L 
with R being the dot (wire) radius and Xso = h/m*a 
the spin-orbit length i 29 ' 42 Then, we can treat the SOI 
within perturbation theory. We assume in the following 
hard-wall boundary conditions for the electrons confined 
in the QDs, namely circular hard-wall boundaries in the 
transverse direction. In the longitudinal direction the 



electron is also confined by a hard-wall type of poten- 
tial, but much stronger than in the transverse direction, 
as stated before. We compute the operator S from Eq. 
© within the first order in SOI, S w (L +L z )- 1 H so , 
which gives explicitly 



(23) 



in the limit of Ez « AE n with AE — Ex — Eq being 
the energy difference between the first excited state |1) 
and the ground state |0). In the above formulas the Li- 
ouvilleans L ,z are defined as L ,zA — [H ,z, A] \IA and 
£ = ^so( — V' x > 0)> ' = B/B. We can obtain an effective 
Hamiltonian up to second-order in SOI and first order 
in Zeeman splitting for the lowest Kramers doublet by 
averaging over the orbital ground state \0t), 



1 



H s „ 7 = -gp B B-<T+(0\[S,H so ]\0) + (0\[S,H e -~ 
1 



-<0|[S, [S,H e ^]]\0)+H,. 



(24) 



The orbital wave-functions have the form (for circular 
hard-wall boundary conditions) 



1 



s/nR J\ m \ + i(k mp R) 



where J\ m \(k mp r) are the Bessel functions of the first 
kind, r is the electron radial coordinate in the trans- 
verse direction, and k mp are the solutions of the equa- 
tion J\ m \(k mp R) = 0. The appropriate energies are given 
by E mp = h 2 k^ np /2m* . Also, we assume that the mag- 
netic field B and the fluctuating electric field E are 
both along the y direction, such that i/e-7 = e£y and 
S = i£ ■ a — (Ez / \so)LQ 1 y a z . After performing the in- 
tegrations, we are left with the following effective Hamil- 
tonian 



1 



H s ^ = -Efa z + + a)a y + H, 



7' 



with 



E 



E z 1-0.25 



( R 



M^eER^. 



(26) 



(27) 



(28) 



We see that there is no second order contribution in SOI 
to the spin-photon interaction, this contribution vanishes 
identically for cylindrical wires in the ground state. We 
mention that within the RWA the Jaynes-Cummings cou- 
pling v becomes v = M?*. 

In the case of two spins present in the cavity, one ob- 
tains the same expression as in Eq. where is 
given by Eq. (|28p . Since our coupling is proportional to 
the bare Zeeman splitting Ez-, we need large magnetic 
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fields in order to obtain a strong coupling. Then, we can 
in principle neglect the Lamb and the ac Stark shifts in 
the expressions for Eg, since they give negligible renor- 
malizations, so that E l z w EfS . However, as can be seen 
from Eq. (|2"?| . the Zeeman splitting can be strongly re- 
duced for large SOI. This feature will turn out to be very 
important in order to have a long-lived qubit (see below). 



V. STRONG TRANSVERSE CONFINEMENT 

In this section we analyze the case shown in Fig. QJ3, 
i.e. when the transverse confinement is much stronger 
than the longitudinal one. As in the previous case, we 
can derive an effective longitudinal Hamiltonian by av- 
eraging the full Hamiltonian H over the transverse or- 
bital ground-state |0t). The effective Hamiltonian H e g = 
Hi = (0t\H\0t) reads 



Hi 



o 
Px 

2m* 



+ V(x)+H z + H l so +H e - 1 +H J , (29) 



with V(x) = (Ot\V (r)\0t) , while all other terms remain 
the same, since they have no action along the x-direction. 
Again, like in the previous case, we disregard the term 
(Ot|(Py + pl)/2m*\0 t ), since it gives a constant shift of 
the levels. 

We now derive the spin-photon interaction from the 
effective Hamiltonian (|29|) . As can be seen from Eq. |3]), 
the SOI contains only one spin-component, along the 
77-direction. In this case and in the absence of an ex- 
ternal magnetic field the SW transformation © can be 
performed exactly, since the SOI appears as an Abelian 
gauge-potentiah 43 i 44 In the presence of an external mag- 
netic field, however, this cannot be done exactly anymore. 
We now apply the transformation © to the Hamiltonian 
Hi so that we obtain Hi = e~ s He s , with the operator 
S corresponding to the zero-field case. This operator S 
reads 



S = -i 



1 ^ . 



with A50 = h/m*r). The effect of this transformation 
can be evaluated exactly and we obtain 



Hi 



Px 

2m* 



+ V{x) + H z {x) + eEx + hiua'<a, (31) 



with 



H z {x) = ig^s^cos {j^ 
fBijCr^ — sin ( -r— ) (e v x B) ■ a 



V Aso 



(32) 



where B v ± is the component of the magnetic field B 
perpendicular to the vector 77, B v is the magnetic field 
component along r?, and e v = r\jr\. We now assume, as 
before, that the Zeeman splitting Ez — S^bB is much 



smaller than the orbital level spacing AEq given by the 
first two term in the above Hamiltonian. Also, we assume 
harmonic confinement potential along the x-direction 
V{x) — iti*ujqX 2 /2 which gives a dot size I — y/hJmFujQ. 
This is usually the case for gate-defined QDs. Then, the 
above condition translates in having Ez <C Hloq. We 
are now in position to derive an effective spin-photon 
Hamiltonian by treating Hzix) within perturbation the- 
ory. We perform a new SW transformation and transform 
the above Hamiltonian into a diagonal one in the basis of 
H a to obtain -ff s _ 7 = (0\e~ s He s |0). We averaged also 
over the orbital ground state |0) to obtain a pure spin- 
photon Hamiltonian. Within lowest order in Ez/hcoo the 
transformation is given by S' = (l-V)L 1 H z (x). After 
inserting the operator S' in the expression for H s _ 1 and 
keeping only the lowest order corrections, we obtain 

-ff s _ 7 = -gfigBeff ■ <T + Al 7 • cr(a^ + a) + froja^a, (33) 

with 

B eff = e-^ Xso ^B nX ■ a + B n a v , (34) 



M a = eVo l J_^ e -m S o) 
d Xso riuiQ 



(er,xl)-a. (35) 



We see that the spin-photon interaction is maximal when 
the magnetic field is perpendicular to rj, like in the per- 
turbative calculation of the previous section. This is ex- 
pected since, as in the previous section, the SOI man- 
ifest itself as an Abelian gauge potential within lowest 
order, although there are two spin-components. For the 
rest of the paper, we assume now a magnetic field per- 
pendicular to 77 so that B^ = 0, B ■ cr^^ = Bcr^ and 
(e,j x /) • a — Ori^.i = ctj- Then, the spin-photon Hamil- 
tonian reads 



# s - 7 = + A4 7 cr £ (a t 



(30) with 



M- 



I I E 



eff 



d Xso hu 



(36) 



(37) 



where E^P = Eze ( ( / A so) 2 is the effective Zeeman split- 
ting. 

We see that the SOI reduces strongly the Zeeman 
splitting for large values of the ratio I /Xso- This over- 
screening of the Zeeman interaction can be understood as 
follows. After performing the transformation (|30|) there 
is no SOI present in the system, but the magnetic field 
in the new 'frame' has an oscillatory behavior, as shown 
in Eq. (f3"2")) . This means that the magnetic field pro- 
cesses around the ^-direction, the speed of precession be- 
ing given by the strength of the SOI measured through 
the SO length Xso- If the bare Zeeman splitting E z is 
much smaller that the orbital level spacing, Ez <C hu>o, 
the electron find itself in the orbital ground state |0) 



7 



given by Hq. Then, if the SOI strength is increased, 
the precession frequency increases also, so that there are 
many precessions of the magnetic field over small dis- 
tances. Since this implies also small changes of the or- 
bital wave- function, this leads to an average reduction of 
the effective Zeeman splitting, as obtained above. 



VI. NUMERICAL ESTIMATES 

We give now some estimates for the coupling v = M v ph 
for QDs in InAs nanowires for both geometries shown 
in Fig. [U In the first case, we assume the dots to 
have a width of 5 — 10 nm (E w w lOmeV) and a ra- 
dius R i=a 50 nm (AEq w 5 meV) . The electron in the 
QD is characterized by m* = 0.023 m e , g ~ 2.5 and 
Xso ~ 100 nm£ We assume also that the ID cavity 
is 2 mm long and 100 nm wide, which implies for the 
fundamental mode Tiuj ~ 0.5 meV and an rms electric 
field E ~ 100 V/m. The Zeeman splitting is assumed 
to be on the same order with the lowest cavity mode, 
i.e. Ef w 0.5 meV (B » 1.75 T). Plugging in all 
the numbers in the formula for v, Eq. (]28[) we obtain 
v w 10~ 5 meV which, in the degenerate case E e f = hui, 
corresponds to a dynamics of the spin-photon system of 
about 60ns (Rabi oscillations between the spin and the 
cavity). In the second case there is more control on the 
orbital level spacing since the dots are obtained in prin- 
ciple by external gating. We now assume a dot radius 
R w 10 nm (E ot ~ 30meV), a dot length I ~ 40 nm 
(Tluj ~ 2 meV) and g « 10. 6 For Ef » 0.5 meV we 
need a magnetic field B i=a 0.45 T. Also, we assume the 
same lengths for the cavity as for the first case so that we 
obtain v « 4 • 10 -4 meV. This gives rise to a dynamics 
of the spin-photon system of about 2 ns in the degen- 
erate limit E e % — Tiuj. We mention that in both cases 
the renormalized Zeeman splitting is quite different from 
the bare one, i.e. E e f — 0.93Ez in the first case and 
E e f = 0.8AEz in the second case. 

For the exchange coupling J between two spins one 
can achieve values as large as J w 10 -6 meV in the limit 
of quite small detunings (A w 10 -4 meV) for the case in 
Fig. [TJ\_, which eventually translates into a time dynam- 
ics of about 500 ns for coherently swapping the two spins. 
In the geometry shown in Fig. [TJ3 the exchange coupling 
J can be much larger, on the order of J sa 4 • 10~ 5 meV 
for detunigs on the order of A sa 4 - 10~ 3 meV, which 
implies a time dynamics of about 20 ns for swapping the 
two spins coherently. 

In order to control the exchange coupling J, one should 
be able in principle to change the Zeeman splitting or the 
orbital level spacing. In InAs QDs the Zeeman splitting 
can be changed very efficiently by changing the dot size 
along the wire direction,— in both cases in Fig. [T] Consid- 
ering the case of two QDs in the cavity, one way to decou- 
ple them is by tuning the g-factors so that Ai = — A2, 
as can be seen from Eq. (fT5|) . However, in the case of 



many QDs inside the cavity this will be rather difficult 
to achieve. 

Another possibility is to change the ^-factors locally so 
that the coupling between the spins reduces due to the 
reduction of the Zeeman splitting Ez- Assuming that 
a reduction of J by one order of magnitude is a good 
measure for the decoupling, one obtains a corresponding 
change in the g-factor of the order of 15% in the first 
geometry shown in Fig. |T|). The rather drastic change 
of g-factor was already experimentally demonstrated for 
InAs QDs by Bjork et al. .6]. They achieved a change in 
the g- factor from \g\ — 3.5 to \g\ — 2.3 when the dot size 
along the nanowire was reduced from 10 nm to 8 nm, i.e. 
a variation of about 30%, which shows to be sufficient for 
our scheme in the geometry shown in Fig. [T] The same 
can be done efficiently for the second geometry, since the 
dots being gate-defined can be modified strongly along 
the wire axis. 

Yet another way to change the exchange coupling J 
is by changing the orbital confining energy AEq. In the 
first geometry v ~ R , and J ~ v (assuming two equal 
spin-photon couplings for simplicity) one obtains a de- 
pendence J ~ R s . Then, by using top gates, for exam- 
ple, one can strongly modify the exchange coupling J by 
a small change of the orbital energy AEo. This can be 
done equally, and maybe more efficiently, for the second 
geometry since, as explained above, the dots can be mod- 
ified easily along the wire axis. The spin-photon coupling 
v ~ Z 4 , which implies then a scaling of the exchange cou- 
pling J ~ I s . 

VII. COHERENT MANIPULATION 

One way to coherently manipulate and to read-out 
(measurement) the qubits is by applying an external driv- 
ing field to the cavity with a varying frequency H ex = 
e(i)(a 1 'e _4 ' Jext +ae 4 " ext ), where e(i) is the amplitude. In 
the dispersive limit (z/j/Aj <C 1) H ex — > H ex + [T,H ex ] 
so that 

H ex ~ etfjate-*"-* + £ ^l a ,+ e — + h.c. (38) 

»=1,2 1 

The control of the i-th qubit can now be realized by tun- 
ing the frequency of the driving field to uj ex — E*% + 
vf I Aj, while this condition is not satisfied for the other 
qubits. This gives rise to an electric-dipole spin reso- 
nance (EDSR) for the i-th qubit, similar to that studied 
by Golovach et ali^ The measurement can be performed 
by tuning the frequency of the driving close to the cavity 
mode so that we can observe peaks in transmission at 
the positions hui + I ^i) a l- If detunings are chosen 

so that all combinations can be distinguished, one can 
measure all the spins from one shot (or at least group of 
spins) 1^ 

A more efficient way to manipulate the spin is to make 
use of the EDSR-scheme proposed in [20( , namely to ap- 
ply an alternating electric field £ (f) to the QD, which, via 
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the electric dipole transitions and the strong SOI, gives 
rise to an effective alternating magnetic field. Briefly, 
if only the dipolar coupling to the alternating electric 
field £(t) is considered, we get H e ^ e i(t) — e£(t)y, with 
the electric field £{t) along y-direction. If the system 
in Fig. is considered, the effective spin-electric field 
coupling within first order in SOI becomes H s - e i = 
(0\[S,H e - e i(t)]\ff) = SB(t)a y , with the fluctuating mag- 
netic field 8B(t) having the form 



SB(t) ~ e£(t)R 



E z R 

AEn A 



*so 



(39) 



For the case shown in Fig. [Tj3 we obtain a similar ex- 
pression for 5B(t), but with the bare Zeeman splitting 
Ez substituted with the effective Zeeman splitting E e ^ 
defined after Eq. ([57)1 , and the radius R substituted with 
the dot length I. The electric field £(t) is assumed to have 
an oscillatory behavior, £ (t) — £ cosuj ac t with u> ac be- 
ing the frequency of the ac electric field. By tuning the 
frequency of the oscillatory electric field oj ac in resonance 
with the qubit splitting E'jf one can achieve arbitrary 
rotations of the spin on the Bloch sphere on time scales 
given by the Rabi frequency ujr = 5B{Q)/h. 2 — We men- 
tion that within lowest order in SOI the induced fluc- 
tuating magnetic field SB(t) is always perpendicular to 
the applied field B and reaches the maximum when the 
applied electric field £{t) points into the same direction 
as This is the reason for choosing the electric field 
along the y-direction. 

We give here also some estimates for the strength of 
the Rabi frequency ojr. For this we assume the same pa- 
rameters as in the previous section and we choose for the 
amplitude of the electric field £q « lOeV/cm. With this 
values we obtain for the strength of the Rabi frequency 
uir w 10 GHz, which gives a time dynamics for the elec- 
tron spin control on the order of lu^ 1 w 0.1ns. This time 
scale must be much shorter than the usual relaxation and 
decoherence times for the spin in the QD. Finding these 
time scales is the subject of the next section. 



VIII. SPIN RELAXATION AND 
DECOHERENCE 

We address now the issue of relaxation and decoher- 
ence of the spin in the cavity. There are two types 
of contributions to the relaxation processes, one aris- 
ing from the finite decay rate of the cavity, k, and the 
other one from the intrinsic relaxation and decoherence 
of the spin, labeled by T± 2 ■ To reach the strong coupling 
regime described here, the losses must be smaller than the 
coupling between the qubits J in the regime of interest 
(y 1 f A > k, Tj ~2 )• Very high-Q factor ID electromagnetic 
cavities were already built (Q = k^ 1 ~ 10 4 — 10 6 ),— so 
that the intrinsic relaxation and decoherence of the qubit 
show up as the limiting factors for reaching the strong 
coupling regime. 



The relaxation and decoherence of the spin-qubit arise 
mainly from the coupling to the bath of phonons and the 
collection of nuclei in the QD. The phonon contribution 
was studied microscopically in great detail for the case of 
gate-defined GaAs QDs in 2DEGs and it was shown that 
for large B-fields, similar to the present case, the main 
contribution to relaxation comes from the deformation 
potential phonons with a decay time T\ ~ 10~ 2 — 10~ 4 s. 42 
As a consequence, a smaller relaxation time is then ex- 
pected for InAs QDs since the SOI is one order of magni- 
tude larger than in GaAs (T\ cx (Xso / R) 2 )- However, 
different from the bulk case, the phonon spectrum in 
nanowires becomes highly non-trivial due to the mixing 
of the branches by the boundaries,— leading to a strong 
modification of the relaxation time. 

In cylindrical nanowires there are three types of acous- 
tic modes: torsional, dilatational and flexural.— All these 
modes couple to the electric charge and, in principle, all 
of them couple also to the spin for a general SOI Hamil- 
tonian. However, as shown later, this is not actually the 
case for the SOI acting in the two configurations in Fig. 
[TJ and only a small part of the entire spectrum gives rise 
to spin relaxation. 

As stated above, within the large Zeeman splitting 
limit considered in this paper, we can take into account 
only the interaction of the electron with the lattice via the 
deformation potential. The clcctron-phonon deformation 
potential interaction is given by H e ^ p h = SoVtt(r,i), 
where Sq is the deformation potential strength and 



u(r,t) 



(40) 



with the displacement field u(k,r) given 

u(k,r) = V$ + (V x e z )$i + (V x V x e z )$ 2 . (41) 

The index k = {q, n, s} quantify the relevant quantum 
numbers, i.e. the wave- vector along the wire, the winding 
number and the radial number, respectively, bk(t) is the 
annihilation operator for phonons, e z is the unit vector 
along the z direction and 



(42) 



with i = 0, 1, 2, n = 0, ±1, ±2 . . . . The functions p n8 (r) 
depend only on the radiu o 45 ' 47 and Xi are normalization 
factors. The effective spin-phonon interaction can be 
found following the same procedure as that used for de- 
riving the spin-photon interaction for both cases in Fig[TJ 



A. Spin-relaxation in strongly-longitudinal 
confined QDs 

We give here the main steps in the derivation of the 
relaxation rate for the case shown in Fig. [TJA Keeping 
only terms up to first order in SOI, we obtain 



H s - P h — (0| [S, H e - p h] 



(43) 
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with S given in Eq. (|23|) and |0) being the orbital ground- 
state. Due to the circular symmetry, the first order in 
SOI term couples only to the n = 1 phonons. The result- 
ing spin-phonon coupling has the form 



H s -ph = -gnB5B y (t)<jy, 



(44) 



with 



C(g,s) « 0.25 / 
./o 



1 dr rJ^fcn^Jofeor)/^^) 
|J 2 (fcn)Ji(fcio)| 



(45) 



, (46) 



where K qs — uj q ^ s / ci with cj 9jS being the eigen-modes of 
the phonon field, c\ the longitudinal speed of sound in 
InAs. The normalization function ^(q, s) is given by 



10 4 

't/2 

H 

io 2 





10' 

1 i J t J u / 

ofzR/ci 

FIG. 2: The relaxation rate T x -1 as a function of the ratio 
u>z R/ci, for both FSCB and CSBC (see text for explanations 
of FSBC and CSBC). Here hci/R ~ 0.6 • 10~ 4 eV (a ~ 4 • 
10 3 m/s and 7? ~ 50 nm) corresponding to a magnetic field 
B ~ 0.2 T, for 3 = 2.5. 



hR 2 



(47) 



where M is the mass of the ions in a unit cell. 

The explicit forms for the u> qs and T{u qjS ) depend 
on the boundary conditions used for the phonon field. 
The two quantities relevant for the boundary conditions 
are the displacement vector u(r) and the stress vector 
t(r) = Te r at r = R, with T being the stress tensor— 
and e r being the unit vector along r. One can now write 
u(r) = Ux and t(r) = Tx with x = {Xo,Xi,X2}, where 
the expressions for the matrices U and T are given in 
Appendix A. There are two limiting cases for the bound- 
aries. The first case is met when there is zero stress at the 
surface, i.e. t(R) = 0,— with ui qjS being the solutions of 
|T(JZ)| = (free surface boundary conditions or FSBC), 
while the second limiting case is met when the surface 
is rigid, i.e. u(R) — 0, with oj q . s being the solutions 
of \U{R)\ = (clamped surface boundary conditions or 
CSBC). The phonon field is normalized according to the 
following relation^ 



2 r27T />R 

^ / d<p / drru* (fc, r, <p) ■ u(k,r,t 

t^R Jo Jo 



h 



2Muj k 



(48) 

From the FSBC or CSBC, together with the normaliza- 
tion of the phonon field, one obtains the spectrum oj q s 
and the normalization function T{q, s). 

We now use the effective spin-phonon Hamiltonian 
with the fluctuating field given in Eq. (|44|) to find the 
spin relaxation and decoherence times, T\ and T2, respec- 
tively. We mention here that the fluctuating magnetic 
field SB y (t) is perpendicular to the external one B such 
that there is no pure dephasing coming from the interac- 
tion of the spin with phonons in lowest order in SOI. In 
fact, as shown previously^ this is valid for any type of 
baths, be it phonons, particle-hole excitations etc. 



In the following we derive the expressions of the T± and 
T2 times resulting from the fluctuating field SB y (t). For 
this we need to compute the bath correlator 



< 8B v (p)6By(t) >, (49) 



where the brackets < ... > means tracing over the phonon 
bath being at thermal equilibrium at temperature T. The 
relaxation time within the Bloch-Redfield approach is 
given in the present particular case (the _B-field along 
x-direction) by (see Ref. 



t; 



MJyy("f) + Jyy(-"f)), (50) 



eff 



Ef jfi. Making use of Eq. (|49| we then 



with u> 

finally obtain for the relaxation rate 
5 



rr 1 



where 



J (0)l 



wfR 
ci 



E 



rp-l 
1 (0)1 



0.05 



duj 



— \ 2 

"0 



S 2 h 



p c R 5 \AE 



C 2 (q,s) 
T{q, s) 



— V 

Aso/ 



(51) 
(52) 



In the above expression 5 — Ez / E e ^ , and the functions 
C(q, s) and F(q, s) are defined in Eqs. (|46I47|) . We men- 
tion that within first order in SOI the decoherence time 
T2 induced by phonons satisfies T2 = 2Ti since, as men- 
tioned before, the fluctuating magnetic field induced by 
phonons SB is perpendicular to the applied one B. In 
Fig. [2] we plot the relaxation time as a function of the ra- 
tio uj e £ R/ci, for R = 50 nm and q = 4 • 10 3 m/s. We see 
that the relaxation rate exhibits peaks as a function of the 
effective Zeeman splitting E'jf . This is due to the finite 
size in the transverse direction which gives rise to phonon 
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branches. Each new peak appears when Eg reaches a 
new energetically higher branch. Note that although the 
relaxation rate seems to diverge when reaching a new 
peak, in reality this does not happen since there are 
many processes which broaden the phonon DOS at these 
special points, like phonon-phonon scattering, phonon- 
substrate scattering, etc. The usual branch splitting is 
on the order of u>p\ = q/ R, which stands for the phonon 
frequency in bulk material with the wave-length equal 
to the dot size R. This frequency ojp\ (or energy, when 
expressed as hu>^ h ) is the parameter which characterizes 
the dominant mechanism for the phonon-induced spin re- 
laxation, which can be due to piezoelectric-potential or 
deformation-potential phonons. In the limit wj? <C 0J ph 
the piezo-phonons give the main contribution to the re- 
laxation rate T-j -1 , while in the opposite case, u>^ 3> ^p\, 
the main contribution to the relaxation rate Tj -1 is given 
by deformation-potential phonons.™ Here we are in nei- 
ther of the two limits, but in the range where Zeeman 
splitting is slightly larger than Tno ph , i.e. u>^ > u> ph . 
However, taking into account only the deformation po- 
tential mechanism should give the right order of magni- 
tude for the relaxation rate. We mention here that the re- 
laxation rate T^ 1 in the low energy limit (uj e £ R/ci < 1) 
is given predominantly by the longitudinal linear in q 
mode (u>i ong (q) = ciq) and the bending mode, square in 
q (<^bend(q) = Bq 2 , with B being a constant which de- 
pends on i?).— 

We see from Fig. [2] that each new phonon branch gives 
a strong enhancement of the relaxation rate , since it 
adds more phonon density of states. However, we see also 
that before the first peak, i.e. before reaching the first 
new branch, there is little spin relaxation (Ti < 10 _3 s) 
for both FSBC and CSBC. This energy scale corre- 
sponds to a Zeeman splitting E'jf ~ 10 _4 eV (E^ fts 
1.2 ■ 10~ 4 eV) for FSBC (CSBC). 

If one tunes the effective Zeeman splitting E% below 
the first peak, the relaxation rate of the qubit becomes 
very small, and the fact that E e ^ and not Ez has to be 
tuned is practically an advantage for reasonably strong 
SOI since we need quite large Ez for having large g oc Ez- 
In the present case E^P /Ez ~ 0.93, and for larger SOI 
this ratio will be even smaller. 



relaxation simplifies considerably. 

The transformation H e ^ p h — > e~ H e - P he , with S 
given in Eq. (|30|). although exact, does not lead to a 
coupling of the spin to the phonon field since both the 
electron-phonon interaction Hamiltonian H e - p h and S 
contain only coordinate x operator, i.e. they commute. 
After this transformation, however, we are left with no 
SOI term, but with the rr-dependent Zeeman coupling 
in Eq. (|32[) . We now perform a second transformation 
H e -ph — > e~ s H e _ p he s with S' given before Eq. (j33j) . 
under the assumption Ez <C Hloq. Then, within first 
order in Ez/nujQ we obtain for the spin-phonon Hamil- 
tonian H s ^ p h the following expression 

H s ^ ph = {0\[S',H e ^ ph }\0), (53) 

where we averaged also over the ground-state |0) of the 
orbital Hamiltonian Hq. The spin-phonon Hamiltonian 
H s - ph reads 

H s -ph = ^9fJ.BSB s .(t)o- i + -gfj, B 6B s (t)(T S) (54) 

with 

hUJ ° qs S)p c UJ q , s /h 

(55) 

and k = {q,s}. The functions M^ ph are given by the 
following expressions 

Mf_ ph (q) = Sinhlnt (J?£) (56) 

Ml ph (q) = 7 - Coshlnt (J^j + Log , (57) 

where 7 = 0.577 is the Euler constant, Log(a;) is the 
natural logarithm, while the special functions Sinhlnt (x) 
and Coshlnt (as) are defined as 

Sinhlnt(x) = f dt^m (58) 
Jo t 



B. Spin relaxation in strongly-transverse confined 

QDs 

We give here a brief description of the phonon-induced 
spin relaxation for the case shown in Fig. [TJ3. We first 
mention that due to the strong confinement in the trans- 
verse direction we can average the electron-phonon inter- 
action over the transverse orbital ground state |0 t ). Since 
for the ground state wave function we have m — (see 
Eq. (|2"5"]l). the only modes which couple to the electron, 
and thus eventually to the spin, are the n = modes 
of the phonon field in Eq. (|4"0")) . Then, the problem of 



Coshlnt(x) = 7 + Log(x) + /" dt ° 0Sh ® . (59) 

We see that, there is both relaxation and pure dephasing 
of the spin due to spin-phonon interaction. However, 
since the deformation-potential phonons is supcrohmic 
(even in ID case for deformation-potential phonons), the 
pure dephasing rate vanishes^ so that we retain in the 
following only the first term in Eq. (|54| . The relaxation 
rate Tj -1 can be found by the same procedure as in the 
previous case and reads 

T^ 1 = ReiJ^iuf) + Jss(-uf)), (60) 
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where the correlation function Jjj is defined in Eq. (|49| 
with y — > x, and to e ^ = Eg /%, as before. The expres- 
sion for the relaxation rate Tf 1 becomes 



Tf 1 = TT 1 



eff 



(0)1 



E 



dq 



q.s 



M 2 s % h (q) 



where 



and 



T, 



(0)1 



2irp c R 2 l 3 \hoJo J 



M* h (q) = M £ s _ ph (q)e 



(61) 



(62) 



(63) 



In order to find now the dependence of the relaxation rate 
Tf 1 on the effective Zeeman splitting ui% , we have to 
find first the phonon eigen-frequencies u) q , s . This can be 
done following the same steps as in the previous section, 
depending which kind of boundary conditions are used, 
i.e. FSBC or CSBC. As mentioned earlier, the average 



R _ 



distance between the branches s is on the order of 
q/R. Then, since R <C I, and also due to the gaussian 
suppression in Eq. (|63|) . it is enough to consider in Eq. 
(|6ip only the lower branch s = 1. If we now assume 
FSBC and the limit qR <C 1, the phonon eigen- frequency 
becomes linear in q, i.e. cj 9j i = uj(q) — cq, with*^ 



'3c?-4cg 



(64) 



The normalization function xa acquires also a simple 
form in this limit, and reads 



R 



Xo 



3c? 



4c^ q V 2Mcq 



(65) 



After inserting in Eq. (|6ip the expressions for ui{q) and 
XOj w e obtain for the relaxation rate Tf 1 (FSBC) the 
final expression 



T, 



rr 1 = 



(0)1 



3cf 



Mf_ ph (u,fl/c) 



(66) 

In Fig. [3] we plot the relaxation rate T, as a func- 
tion of the dimensionless parameter uj^I/c for different 
SOI strengths measured through the ratio l/Xso- We 
assumed here R = 10 nm and I = 50 nm, which gives 
kc/l ee TvJ vh = 0.05 meV and Tici/R = huj* h = 0.25 meV. 
We see in Fig. [3] that the relaxation rate T^ 1 is quite 




3 4 

(£>fl/c 

FIG. 3: The relaxation rate T^ 1 as a function of the ratio 
lj^I/c for three different ratios l/Xso and with FSBC (see 
text). 



branch since above it we find a substantial increase for 
the relaxation rate. Since this next phonon branch 
lies somewhere around 2%ojp h « 0.5 meV, the condi- 
tion for strong spin-phonon coupling and weak relax- 
ation becomes Tiuj l ph <C < 2fwj ph . In this regime we 
are actually satisfying also the necessary condition that 
Ez /hu>o <C 1 , since for I = 50 nm we have Tiloq = 1.3 meV. 
We mention that for CSBC the phonon spectrum is 
gaped, and, in consequence, there is no phonon-induced 
relaxation of the spin for Zeeman splittings Eg smaller 
than this gap A p h- This energy (gap) is on the order of 
Ap/j ~ 2hWph — 0.5 meV. Note the non-monotonic be- 
havior of the relaxation rate as a function of the effective 
Zeeman splitting (see Fig. [3j. This non-monotonicity 
has the same origin as in GaAs QDs^ and comes from 
the fact that for increasing Zeeman splitting the wave- 
length of the phonon decreases, and when this becomes 
less than the dot length the phonons decouple from the 
electron (i.e. the electron-phonon coupling averages to 
zero). A similar non-monotonic effect has been recently 
observed in GaAs double QDs. 51 

The spin decoherence time due to phonon processes 
is given by T2 = 2T\ such that the main source for de- 
coherence comes from the hyperfine interaction between 
the electron and the surrounding nuclei. This time scale, 
T 2 * , is given by^^ 



To 



N 

X 



(67) 



large (Tf 1 



10 5 



where N is the number of nuclei in the sample and A 
is the hyperfine constant. We see that the larger the 
number of nuclei, i.e. the bigger the dot, the longer is the 
pure decoherence time for the electron. In a typical 
GaAs QDs (R = 30 nm and I — 5nm) this time scale 



10 7 s 1 ) for u'f /uj l h — 1 - 5, i.e. is on the order of T 2 * — 10~ 8 Sr2£2 which corresponds 



when these energies are comparable. However, there is 



need for a large effective Zeeman splitting E'jf ^> Hlu 
to achieve a large spin-photon coupling M. 1 . At the 



pit 



also to our two cases described in the paper. However, 
again like in GaAs, we expect that coherently driving 
the qubit will prolong the T 2 * time up to 10~ 6 s and with 



same time, one should stay still below the next phonon echo up to 10 sJ* Moreover, like in GaAs QDs, one 
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can make use of state narrowing procedures ; > 55 which 
should lead to a further substantial enhancement of T 2 * 
due to nuclear spins, and possibly reach the SOI induced 
limit of 10 _1 — 10 -4 s calculated above for large magnetic 
field strengths. 



IX. CONCLUSIONS 

We have proposed and studied an efficient way to im- 
plement spin qubits localized in InAs nanowires cou- 
pled to a ID electromagnetic transmission line (cavity) 
via SOI. We have analyzed two experimentally achiev- 
able configurations of the system. In the first case 
the electronic confinement is much stronger along the 
nanowire axis than in the transverse direction (large- 
radius nanowires), while the other case corresponds to 
the opposite limit (small-radius nanowires). We have 
found a reasonably strong coupling between the spin and 
the cavity modes due to both strong vacuum fluctua- 
tions in the cavity and strong SOI in InAs. We also have 
shown that this spin-photon coupling can allow for cou- 
pling between two (or several) distant spins, depending 
on the detuning of the Zeeman splittings E e ^ from the 
cavity mode hu. The SOI- induced exchange coupling J 
between two spins can be controlled by electrical fields 
only, e.g. by changing the g- factor and/or orbital level 
spacing. Also, single-spin rotations can be performed ef- 
ficiently by electric fields only, through the EDSR mech- 
anism. In principle, the price one has to pay in strong 
SOI materials is strong coupling to the charge environ- 
ment which then relaxes and decoheres the spin. How- 
ever, we have studied the relaxation of the spin due to 
the lattice vibrations in the InAs nanowires for both con- 
figurations, and shown that the time scale for the spin- 
decay is on the order of milliseconds for strong magnetic 
fields (B ~ 0.5 —IT), much larger than the times asso- 
ciated with the spin-photon dynamics, which takes place 
on times scales on the order of 10~ 8 — 10~ 7 s. This fact 
is due to the quasi ID structure of the system where the 
phonon spectrum shows discrete branches, very different 
from the bulk limit. 

We stress here also that the coupling of the quantized 
modes of the transmission line to the spin degree of free- 
dom via SOI is not restricted to QDs in semiconduc- 
tor nanostructures. In principle, this coupling should 
be possible in other spin-orbit coupled systems too, 
like nitrogen- vacancy centers (NV-centers)) 56 i 57 molecu- 
lar magnetS ) 58 ' 59 i 60 magnetic nanorings^i etc. In these 
systems there is usually a strong zero-field splitting (ZFS) 
of the lowest spin-multiplet attributed to SOI or to 
dipole-dipole interaction. This would allow for an effi- 
cient coupling of the electric fields, quantum or classi- 
cal, to the spin degree of freedom and finally providing 
a mechanism for an all-electrical implementation of spin- 
based quantum information processing. 

As a final remark, we mention that the present scheme 
can be also used to form hybrid structures where spin- 



qubits are integrated together with other types of qubits 
in the same ID transmission line. For example, one can 
envision a setup where a spin-qubit is coupled via the 
cavity modes to superconducting qubit as the one stud- 
ied in Ref. [35j so that one can transfer arbitrary states 
between the two qubit-systems. 
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APPENDIX A: DISPLACEMENT AND STRESS 
TENSOR FOR CYLINDRICAL NANOWIRES 

In this Appendix we give explicit formulas for the dis- 
placement u(r) and stress t{r) vectors, respectively. We 
can write the displacement vector u(r) — (u r ,u^,u z ) 
from Eq. (|41[) in components 

Ui (r,t) =J2 U v(r)Xje l ^ z -^\ i = r,0,z, (Al) 

3 

with Xi — (Xfb Xi) X2) and the matrix U(r) having the 
form 

/ ■§pfon(r) iffin(r) iq£fo n (r) \ 
U{r) = i^f Qn (r) -&/ ln (r) -**/ 2 „(r) . (A2) 
V iqfon{r) fcfAn(r) / 

The other relevant quantity for the elastic problem is 
the stress tensor T4£ In order to obtain T, we first 
have to find the strain tensor S as a function of dis- 
placement u(r). The independent components of the 
strain tensor coordinates have expressions (in cylindri- 
cal coordinates)^ of the form 

dUr_ 

or 

1 { duj, \ 
S ** = r{-W +Ur ) 

c _ du z 

OZ 

1 fdu r o d fu r \\ 
br4> ~ Yr\~o^ +r d~rV7)) 

s - 1 9uz 1 ^± 

z ^ r dcj) dz 

*--K£ + £)- (A3) 

The stress tensor, T , which quantifies the surface forces, 
is related to the strain tensor S by the elastic modulus 
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constants.™ Since we are interested in the boundary con- 
ditions at the surface of the cylinder, the relevant part 
of the stress tensor is given by the stress vector t = Te r , 
with e r being the unit vector along the radius. We write 
here only these relevant parts of the stress tensor T as a 
function of the strain tensor components 



T, 



r<j> 



p(cf - 2c 2 t )(S r 
2pc t S rt p 
2pc t S rz , 



S Z z) + 2pcfS r 



(A4) 



and t = (T rr , T r( p, T rz ). We write now the relevant stress 
vector t, which is given explicitly by the following relation 



/ T rr 




T r( f, 


H 


\T rz 




stress 


matrix 



+{#-24)1 (c?-2e?U£ (c 2 



„21J. 

-t r dd> 

C t dz 



as we did for the displacement, namely ti(r) 



-*t ( Or r ) 









«2 d_ 

-t dr 




(A5) 



explicit form 



with the matrix T having the 



T(r) 



g r 2 {cf 



2c 2 ; 



77 h 



2 «c 2 g|:/o„ 



2mc 2 



JLfii*.') 

dr \ r J 

(2^ + fc 2 )/ ln -2qnc$jL(ff-) 



\ 



(A6) 
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